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' , Abstract: We investigate the properties of hot and/or dense matter in QCD-hke theories 

QJ \ with quarks in a (pseudo)real representation of the gauge group using the Nambu-Jona- 

'"^ ' Lasinio model. The gauge dynamics is modeled using a simple lattice spin model with 

nearest-neighbor interactions. We first keep our discussion as general as possible, and only 
^ ■ later focus on theories with adjoint quarks of two or three colors. Calculating the phase 

OO , diagram in the plane of temperature and quark chemical potential, it is qualitatively con- 

Q"^ , firmed that the critical temperature of the chiral phase transition is much higher than the 

^N I deconfinement transition temperature. At a chemical potential equal to half of the diquark 

^O [ mass in the vacuum, a diquark Bose-Einstein condensation (EEC) phase transition occurs. 

Q ' In the two-color case, a Ginzburg-Landau expansion is used to study the tetracritical be- 

havior around the intersection point of the deconfinement and EEC transition lines, which 
are both of second order. We obtain a compact expression for the expectation value of the 
^ . Polyakov loop in an arbitrary representation of the gauge group (for any number of col- 

H . ors), which allows us to study Casimir scaling at both nonzero temperature and chemical 

potential. 
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1. Introduction 

Quantum chromodynamics (QCD) is accepted as the theory for strongly interacting matter. 
However, due to the strong coupling, perturbative treatments fail at an energy scale of the 
order of Aqcd , resulting in the fact that the structure of the QCD vacuum is still far from 
being well understood. The same is true of the QCD phase diagram where, due to the 
formidable sign problem, standard Monte-Carlo techniques based on importance sampling 
cannot be used at nonzero quark chemical potential. 

In order to get deeper insight into the behavior of dense quark matter, several QCD- 
like theories have been proposed, including QCD with an imaginary chemical potential [l|, 
at nonzero isospin density Q, two-color QCD with fundamental quarks [^, and any-color 
QCD with adjoint quarks 0]. All these theories share the pleasing feature that they are 



free of the sign problem, which makes them possible to be simulated on the lattice from 
first principles. 

In the present paper we focus on two wide classes of QCD-like theories: those with 
nonzero baryon chemical potential and quarks in a real or pseudoreal representation of 
the gauge group. For the sake of brevity, they will be henceforth referred to as type-I 
and type-II QCD-like theories, respectively (see, e.g., Refs. [^, |6| for applications of these 
theories in another, electroweak, sector of the standard model of elementary particles). The 
typical examples are QCD with adjoint quarks of two (aQC2D) or three (aQCD) colors 
for type I, and two-color QCD with fundamental quarks (QC2D) for type II. While the 
effective Nambu-Jona-Lasinio (NJL) model description for the latter was worked out in 



Refs. ||7|, g, ^, |T^, [ll[], the model Lagrangian for type-I theories is, as far as we know for 
the first time, constructed here.^ Due to the (pseudo)reality of the quark representation 
of the gauge group, all these theories have several remarkable differences as compared to 
ordinary QCD, besides the absence of the sign problem itself. 

First, with A'^^ massless quark flavors, the global flavor symmetry is SU(2A'^^) rather 
than the usual chiral group S\J{Nf)i^ x S\J{Nf)^ x U(1)b. The reason is that the charge- 
conjugated quark field (V'r) (charge conjugation being defined as tp = C'ip ) which 
is a left-handed spinor transforms in the same way as the left-handed quark ^i^ under 
both color and Lorentz transformations, so it is allowed to transform them into each other 
while keeping the color symmetry intact. This means that the multiplets of states in 
the spectrum will contain modes of different baryon number. In particular, apart from the 
pions the Nambu-Goldstone (NG) bosons of the spontaneously broken fiavor symmetry will 
also include diquarks. These light diquarks are colorless bosons carrying baryon charge, 
and hence at low temperature and sufficiently high chemical potential, they will undergo 
Bose-Einstein condensation (BEC). 

Second, in case of quarks in a real (such as the adjoint) representation, the Zn^ center 
symmetry remains intact in the presence of dynamical quarks. This leads to a well-defined 
deconfinement phase transition, accompanied by spontaneous center symmetry breaking. 



instead of a crossover as in QCD |13|. The associated order parameter is the expectation 
value of the Polyakov loop. For the two- and three-color cases investigated in this paper, 
the deconfinement transition is of second and of first order, respectively. 

Since the BEC and deconfinement phase transitions are both well-defined, being asso- 
ciated with exact symmetries even in the presence of massive dynamical quarks, aQC2D 
exhibits a rather unusual critical behavior in the vicinity of the tetracritical point where 



the two second-order transition lines cross each other |14]. In (three-color) aQCD the de- 
confinement transition is of first order. As a consequence the second-order BEC critical line 
is interrupted around the deconfinement transition, meeting the deconfinement line at two 
tricritical points. This general expectation is confirmed by our explicit model calculation. 
To model the gauge sector, we use a simple lattice spin model with nearest-neighbor 
interaction, inspired by the strong-coupling expansion [^, ^, 0, |l8|, 19, 20|. This is 
then coupled to continuum quarks in a fashion similar to the Polyakov-loop NJL (PNJL) 

^A Polyakov loop NJL-type model for adjoint fermions was already worked out in Ref. |12|, but the 
quark sector of their Lagrangian does not have the S\J{2Nf) symmetry of the underlying gauge theory. 



model [^, 21, 22, ^, 24, 25 1. As an effective model, the NJL model successfully describes 
chiral symmetry breaking and pairing of quarks. But since it contains no dynamical gluons, 
the confinement feature is missing. In order to account for confinement at least in a heuristic 
way, one adds to the thermodynamic potential an effective potential for the Polyakov loop, 
adjusted to reproduce the thermodynamic observables in the pure gauge theory. The 
Polyakov loop is represented by a constant temporal background gluon field which in turn 
couples to the quarks [^. The parameters of the PNJL model are fixed separately in 
the pure gauge and chiral quark sectors. The successful qualitative reproduction of the 
coincidence of the deconfinement and chiral restoration temperatures, T^ and T^, in QCD 
is then one of the great virtues of the PNJL model. On the other hand, aQCD is very 
different. First, T^ <^ T^, resulting in a broad range of temperatures exhibiting deconfined. 



but still chirally broken matter |g^, ^, 28 1 (see also Refs. [^, 29 1 for related theories with 
periodic boundary conditions for quarks). Second, Td does not change much compared 
to the pure gauge theory when quarks are coupled in, because adjoint quarks carry zero 
center charge. We confirm these features in our results. 

The paper is organized as follows. In Sec. ^ we introduce our model, working out 
separately the actions in the gauge and quark sectors. The gauge part is well known 
from literature, and we therefore just elaborate on the Weiss mean-field approximation 
used in this paper. In the quark part we deal with the task to construct an interaction 
Lagrangian with S\J{2Nf) fiavor symmetry. While this was previously achieved for QC2D 
and actually applies equally well to all type-II theories, here we construct analogously a 
model Lagrangian for type-I theories. Section ^ is devoted to two-color QCD. We study 
the phase diagram of aQC2D and derive the Ginzburg-Landau (GL) theory that governs 
the behavior of the system near the tetracritical point. We find a simple closed analytic 
expression for the expectation values of the Polyakov loop in all representations, valid in 
pure gauge theory as well as with dynamical quarks in an arbitrary representation. In 
Sec. 1^ we show analogous results for aQCD. Finally, in Sec. ^ we summarize and conclude. 
Some technical details are deferred to the appendices. Throughout the paper, we use the 
natural units in which the Planck's and Boltzmann's constants as well as the speed of light 
are equal to one, and the timelike metric in the Minkowski space. 

2. Model setup 

In this section we set up the model that we later on use for numerical computations. In the 
gauge sector we employ a simple lattice-inspired model, which can in principle be used for 
any number of colors. The quark NJL Lagrangian derived afterwards is, as already stressed, 
applicable to all QCD-like theories with quarks in a real representation. This is natural: 
the Lagrangian is based almost exclusively on the fiavor symmetry and is therefore valid 
for an arbitrary number of colors. The numerical values of the parameters in the model 
will be fixed in the following sections when we come to the discussion of concrete results. 

2.1 Gauge sector 

Our starting point for the pure gauge sector is an effective theory for the Polyakov loop 



inspired by the lattice strong-coupling expansion. We closely follow the notation and line 
of argument of Ref . pO| . The action of the model is given by 



^g[L] = -Nle-''/^ Yl ^^i^yU^ + y) ' (2-1) 



x,y 

where x are the lattice sites and y are the neighboring sites. (We use boldface to indicate 
spatial vectors.) The only adjustable parameter a is related to the string tension and can 
be extracted from numerical simulations of the full (pure) gauge theory. Furthermore, 
i-F{x) = j^ Tr Lp{x) is the traced Polyakov loop in the fundamental representation; in the 
full gauge theory, the Polyakov loop in a given representation 7^ is defined as 



Lti{x) = Vexp 



l/T 

dTAl{x,T)Tan 




(2.2) 



where Ta-jz are the gauge generators in this representation. 

In the so-called Polyakov gauge where temporal gluons have constant values, this 
simplifies to 

Ln{x) = exp [iAl{x)Tan/T] . (2.3) 

Moreover, only the components of A^ corresponding to generators that form the Cartan 
subalgebra of the gauge group are nonzero. Let these components be 6iT. (There are Nc — 1 
independent ones; the conventional factor T makes the variables Oi dimensionless.) Each 
representation of the gauge group is characterized by a set of weights, Wia, that represent 
the eigenvalues of the generators of the Cartan subalgebra in this representation; the index 
a labels the different eigenvectors of the Cartan subalgebra. The traced Polyakov loop in 
representation TZ then reads 

£^(a;) = —1— y e^^'(^)'"- . (2.4) 

a 

In the fundamental representation, the Polyakov loop (in the Polyakov gauge) is usually 
represented as diag(e'^i, . . . , e^^^<:-i, e~'(^i^ ^^iVc-i)). This corresponds to the choice of 
the Nc weights of the fundamental representation as Wia = Jjq for q = 1, . . . , A''^ — 1, and 
WiNc = —1 for all i. Equivalently, it can be written by defining 9^^^ = —{6i + • • • + 6nc-i) 
up to an integer multiple of In. 

In the Weiss mean-field approximation, the nearest-neighbor interaction is linearized 
and the action (|2.l|) is replaced with the action o?%f(a,/3), depending on two mean fields 
a,/3,2 

=^mf(a,/3) = -NcY^[a Re ipix) + i^lm If (x)] . (2.5) 

X 

The dynamical variables of the model ( ^.ij ) are the (untraced) Polyakov loops L{x) and its 
partition function is therefore obtained as 5^ = exp{—Qg/T) = j Y\^ dL{x) exp(— ^^^^[L]), 
where dL is the group-invariant (Haar) measure on the SU(A'^c) gauge group. For the sake 



^Here, we adhere to the notation introduced in Ref. pd]. Let us therefore just stress that the symbol /? 
is not to be confused with the inverse temperature. 



of future reference, let us add that in terms of the phases 9i, the Haar measure can be 
written as 

Nc-l Nc 

dL= l[ d^ij||e^^»-e^^^f , (2.6) 

i=l i<j 

The integration over the variables 9i is performed over the range [0,27r]. 

The thermodynamic potential can now be rewritten by subtracting and adding the 
mean- field action, resulting in the expression 



n, 



log(e-(-^«-^-f)>^f-log /'j]dL(a;)e-^-f. (2.7) 



Here and in the following, {■)mi is the average with respect to the distribution defined by 
the mean- field action. For a given (not necessarily local) function 0[L] of the Polyakov 
loop, it reads 



YldL{x) 0[L]e' 



-y^ 



(0)mf= J . (2. 



ndL( 



X ) e '^™f 



Note that when the function O is local and does not depend explicitly on the coordinate, 
the product over lattice sites can be dropped. 



Equation (2/7) is still exact; no approximation has been made so far. By the same 
token, the thermodynamic potential Vlg is independent of the arbitrary variables a, /3. In 
the Weiss mean-field approximation, one replaces {|e~( s"'^™") . with e~''^s-'^mt)mf po[| . 
The mean fields are then determined selfconsistently from the stationarity condition. In 
fact, as long as /3 = so that the averaging is done with a real mean-field action, one can 
use Jensen's inequality^ to show that this approximation provides a strict upper bound for 
the exact free energy. Its optimum estimate is then obtained by minimizing with respect 
to a. 

The final formula for the Weiss mean field gauge thermodynamic potential reads 

_^ = - 2(d - l)7V,2e-'^/^(£F>mf(^F>mf + ^ [(« + /3)(^F>mf + (o - /3)(£^)n,f] - 

^^ , 2 (2.9) 

— log / dLe^'=*^"^'^^^+^'^^™^^^ . 

Here a^ denotes the lattice spacing and the factor a^/y is just the inverse of the number 
of lattice sites; d stands for the dimensionality of spacetime so that 2{d — 1) is the number 
of nearest neighbors on a cubic lattice. 



^Jensen's inequality (see e.g. Ref. bol]) states rather generally that for any real convex function /, 
f{{x}) < (/(a:)), where the averaging involves either a (weighted) arithmetic mean in the discrete version 
of the inequality, or an integral average over a given probability distribution in the continuous version. 



2.2 Quark sector 

The Lagrangian of the quark sector cannot be derived from the underlying gauge theory 
directly. However, it is strongly constrained by the requirement that it inherits all the 
symmetries of the QCD-like theory. As already stressed above, in theories with A'^^ massless 
quark flavors in a (pseudo)real representation of the gauge group, the usual chiral symmetry 
is promoted to SU(2Aj). In order to see how this comes about, let us start from the 
Lagrangian of the gauge theory, including a common mass mo for all quark flavors, 

^QCD-like = V^i^V' - mo^i^ , (2.10) 

where !^^tp = {d^ — i^TaApV is the gauge-covariant derivative. Indices are suppressed so 
that this formula holds for quarks in any representation of the gauge group. 

The fact that the quark representation is (pseudo)real means that there is a unitary 
matrix ^ such that ^^ has the same transformation properties under the gauge group 
as -0. It is then advantageous to trade the Dirac spinor, consisting of the left- and right- 
handed components, for the purely left-handed Nambu spinor, 

A crucial fact known from the theory of Lie algebras is that ^ is either symmetric or 
antisymmetric according to whether the quark representation is real or pseudoreal pT| . 
Writing collectively ^'^ = ±^, we can introduce the charge-conjugated Nambu spinor. 

The Dirac conjugate of both *& and ^® is defined naturally by conjugating the individual 



components. The Lagrangian ( 2.10 ) then becomes, in the Nambu formalism, 



^QCD-like = ^i^* 



1 -^/ 1\ ^, 

2 " \±1 Oi 



(2.13) 



First of all, we can see that in the chiral limit, the Lagrangian of a QCD-like theory indeed 
has an SU(2Aj) symmetry. Note that baryon number is already incorporated in this 
simple group, for it is represented by the block matrix 2diag(l, —1) in Nambu space. The 
change of the overall phase of the Nambu spinor corresponds to the axial U(1)a symmetry 
which is broken at the quantum level by instanton effects. Since the mass term has the 
same structure as the chiral condensate, we can also immediately infer that for type-I 
(type-II) theories the order parameter for flavor symmetry breaking transforms as a(n) 
(anti)symmetric rank-two tensor of SU(2Aj). Therefore, the two classes of theories have 
different symmetry-breaking patterns and subsequently also different low-energy spectra. 
The symmetry-breaking patterns in the vacuum are SU(2Aj) — )■ S0{2Nf) and SU(2Aj) — )• 
Sp(2Aj) for type I and type II, respectively Q. 



The task to construct an NJL-type interaction compatible with the S\J{2Nf) symmetry 
is most easily accomplished using the Nambu notation ( 2.11| ). It is useful to stress right at 



the outset that as long as only color-singlet channels are considered, each of the Lagrangians 
to be constructed below applies to the whole class of QCD-like theories (type-I or type-II), 
regardless of the detailed structure of the gauge group or the quark representation. In fact, 
NJL Lagrangians for type-II theories with two quark flavors were already constructed in 
Ref. [|l^]. Here we follow the same line of argument with the necessary modifications for 
the type-I case. 

One property that further distinguishes the type-I and type-II theories is the severity 
of the sign problem. While we remarked before that all QCD-like theories considered in 
this paper are free from the sign problem, one should be a bit careful with the type-II 
theories. There, the determinant of the Dirac operator is in general real, but needs not 
be positive. In order that there be no sign problem, one therefore has to consider an even 
number of flavors. On the other hand, type-I theories have no sign problem for any number 
of flavors |32|. As a warm-up exercise, we thus start with the simplest case of one flavor. 



In the following, the Pauli matrices o"o,i,2,3 = {^,(^1,(^2,(^3} are used to denote the 
block matrices in Nambu space, and to,i,2,3 = {^,ti,T2,T3} are used to denote the fla- 
vor generators for Nf = 2. The symmetric rank- two tensor representation of the flavor 
SU(2) ~ S0(3) group is real and three-dimensional. Using the basis of symmetric unimod- 
ular unitary matrices as S = {l,i(Ti,io"3}, we can immediately construct two four-fermion 
interaction terms, 

rr o 1 (2-14) 

Lz / _ .^^ _\9 _.r,— ,0 ,— .0 . . .^ . o . . .^ .0-1 V / 



■^if,su(2) 2 L 



While the former preserves the axial U(1)a, the latter breaks it explicitly. It is easy to 
verify that ^if,su(2) is the 't Hooft determinant term, i.e. 

^if,su{2) = 2G{detWf-^j+R.c.) . (2.15) 

For two flavors, the ten basis matrices of the symmetric rank-two tensor representation 
of the flavor SU(4) are chosen as the symmetric Kronecker products of a and r, i.e. 

■^ — l<7sym '<>' '^symi '^'antisym '<> ''"antisymj • K^'^^) 

Since the 10-dimensional representation of SU(4) is complex, only one of the above two 
possibilities to construct an invariant interaction term remains, 

-^2f,U{4) = G\W^^\' 

= g[(V5V)' + (^i75rV)' + (^iTsV')' + ii^rijf + ^ \Wrsi^\^ + Yl I^TsTs^P 

s s 

(2.17) 

which preserves U(1)a automatically. (Here ts denotes the set of symmetric Pauli matrices, 
Ts = {l,ri,r3}.) A U(1)a breaking interaction can again be introduced by the 't Hooft 
determinant term, but such a term will be an eight-fermion contact interaction which we 
do not consider in our model. 



2.3 Mean-field approximation 

We will employ the usual mean-field approximation, introducing the collective bosonic 
fields via the Hubbard-Stratonovich transformation and subsequently replacing them with 
their vacuum expectation values. To that end, however, one first needs to guess which 
condensates (order parameters) will appear in the phase diagram. The case of type-II 
theories with two quark flavors was worked out in Ref. |jll|: as long as just the baryon 
chemical potential is considered, one only needs the chiral condensate, a = —2G{il)ip), and 
the scalar diquark condensate, A = 2iG{il! C^^^T2'4')- Since the diquark wave function 
is antisymmetric in color as well as spin indices, it must, by means of the Pauli principle, 
also be antisymmetric with respect to flavor. The (spin-zero) diquark in type-II theories 
therefore mixes quarks of different flavors. Consequently, in the presence of an isospin 
chemical potential the diquark pairing feels stress and eventually diminishes via a first-order 
phase transition, with a narrow window of chemical potentials featuring inhomogeneous 



pairing |11 



In type-I theories the scalar order parameters are symmetric in color and antisymmetric 
in spin indices, hence they must be symmetric in flavor. This is in accordance with the 
fact that for two flavors, there are altogether nine NG bosons of the SU(4)/SO(4) coset, 
the isospin triplet of pions and the isospin triplet of (complex) diquarks. At zero isospin 
chemical potential, the isospin multiplets are strictly degenerate. In particular all uu, 
dd, and ud + du diquarks can condense when the baryon chemical potential exceeds their 
common mass. However, for arbitrarily small isospin chemical potential, the diquarks 
formed from quarks of the same flavor will be favored. Such single-flavor condensates do 
not feel stress at nonzero chemical potential, and the phase diagram of type-I theories will 
therefore not contain inhomogeneous phases, as observed in Ref. |34]. 



With the above argument in mind, we restrict our attention to single-flavor conden- 
sates. The fact that the two- flavor four-fermion interaction ( |2.17 ) automatically preserves 



U(1)a means that the condensates differing just by opposite parity will be degenerate. 
However, we know from the Vafa-Witten theorem that in the vacuum parity is not spon- 
taneously broken |35]. The degeneracy will be eventually lifted by instanton effects, man- 



ifested in the eight-quark 't Hooft interaction term. Within the present model, we will 
simply ignore the negative-parity channels. 

As long as we only deal with one-flavor condensates, we can write down the contribution 
to the thermodynamic potential from a single quark flavor. This equals the thermodynamic 
potential of free fermionic quasiparticles. In presence of a pairing gap A, their dispersion 
relation reads Ef. = \/(^£)^ + A^, where ^| = efc + e/x, e = ±, and e^ = Vk"^ + M"^; 
M = mo + cr is the constituent quark mass and /i the quark chemical potential. The 



gauge and quark sectors are coupled in the PNJL spirit |21|. In the Polyakov gauge the 
temporal component of the gauge field is constant. The individual quark color states in a 
given representation will then have, in the presence of the background gauge field, effective 
chemical potentials iT'^-OiWia- Since the quasiparticle spectrum discussed above is the 
same for all color states in the representation (this is because all condensates are color 



singlets!), the thermodynamic potential of one quark flavor will simply be 



VNf AG ^^ J {2ttY 

1/2 (2.18) 

xiEldiian + 2T\oglT\\l + 2cos{eiWio,)e~^'k/T j^(,-^Ei/T \ 1 

a 

The power of 1/2 in the second line compensates the doubling of the number of degrees of 
freedom in the Nambu formalism. 

The group average must be performed once we couple the quarks to the Polyakov 
loop. Note that we do not average the full quark thermodynamic potential, but only 
the argument of the logarithm. This replacement was introduced in Eq. (13) of Ref. [|20| 
as a convenient approximation to (ilg)mf. However, in Appendix O we present a heuristic 
argument showing that the prescription ( p. 18 ) is actually superior to the full average (riq)jnf. 



While with fundamental quarks considered in Ref. |20[ | the numerical difference between 
the two ways of evaluating the quark sector thermodynamic potential is negligible, we point 
out that with adjoint quarks, taking the average (f^g)mf would lead to unphysical artifacts 



which are not present in Eq. ( 2.1 



2.4 Parameter fixing in the quark sector 

The NJL part of the model has three adjustable parameters: the coupling G, the current 
quark mass mo, and the ultraviolet cutoff that regulates divergent integrals. (Within this 
paper, we will use the three-momentum regularization scheme.) These need to be fixed by 
fitting to three selected observables. A conventional, and convenient, choice are the chiral 
condensate, pion mass, and decay constant in the vacuum. While the pion mass is more or 
less a free parameter that can be easily modified in lattice simulations by tuning the quark 
mass, the remaining two parameters depend on the single physical scale of the underlying 
theory, and cannot therefore be adjusted at will. 

In three-color QCD with fundamental quarks, one can directly use experimental observ- 
ables. In QC2D, the input parameters were determined in Ref. [|lO|| from their three-color 
counterparts by A'^c-rescaling. Unfortunately, we are not aware of suitable lattice data that 
would allow us to fix the parameters directly in the case of aQCD and aQC2D. We therefore 
use the following indirect argument. Suppose that we have a theory with both fundamental 
and adjoint quarks. Gauge invariance can then only be maintained when the coupling of 
quarks to gluons is the same in both representations. Since the effective meson-channel 
Lagrangians of the NJL type can be derived from a one-gluon-exchange-inspired interac- 
tion, this allows us to fix the ratio of the effective couplings in the fundamental and adjoint 
quarks sectors. 

Concretely, assume the current-current interaction 

^int = -gii^rTani^f . (2.19) 

The coupling g can be directly related to the microscopic QCD coupling and the screening 
mass of the gluon in the one-gluon-exchange approximation. We therefore assume that it is 



a [MeV] 


6^/3 [MeV] 


A [MeV] 


G [GeV-2] 


mo [MeV] 


670.9 


269.2 


657 


25.71 


5.4 



Table 1: Model parameters for two-color QCD with adjoint quarks. 

the same for fundamental and adjoint quarks. Performing the Fierz transformation to the 
meson channel yields the effective NJL coupling Gp = gi^c ~ ^)/i'^^c^f) f°^ fundamental 
quarks [^]. For adjoint quarks we analogously obtain Ga = gNc/[i^c~'^)^f]- This results 

in the ratio 

Ga 2NJ} 

-^ = . (2.20) 

For the reader's convenience, the derivation of this relation is sketched in Appendix W. In 
the following sections, we will use it to infer the value of the coupling for adjoint quarks 
from that for the fundamental ones. We will not refer to the original current -current 
interaction anymore. 

Eq. ( p. 20] ) would at first glance suggest that the coupling for adjoint quarks is weaker 
than for the fundamental ones (with the exception Nc = 2). One may then wonder why 
the chiral restoration temperature is much higher for adjoint quarks. The reason for this 
is that in the gap equation, the coupling is multiplied by the number of quark degrees 
of freedom coming from the quark loop. The effective coupling ratio for adjoint versus 
fundamental quarks therefore is 2N^ /{N^ — 1) which is always larger than two. 

3. Two colors 

For two colors, the group integration is easily done and it is possible to find closed analytic 
expressions for all general formulas derived above. First, there is just one independent phase 
6, associated with the only diagonal generator of the SU(2) gauge group. The (2j + 1)- 
dimensional spin-j representation then has weights —2j9, • • • , +2j9, and one immediately 
obtains 

^ 2j + l sin6' ^ ' 

The Haar measure (|2.6| ) reduces to dL = ^ sin^ 6 dO, normalized so that the group volume 
is unity. 

Since all traced Polyakov loops of SU(2) are real, we need just one mean field a in 
Eq. ( |2.5| ). Using the definition of the modified Bessel function of integer order, 

/„(x) = i /" d^e^^^^cosn^, (3.2) 

and the recurrence relation /„_i(x) — In+i{x) = ^In{x), one derives the expectation value 
of the Polyakov loops [37], 

hj+ii2a) 



tj/mf 



Ii(2a) 



(3.3) 



10 



The gauge part of the thermodynamic potential (2^) in turn becomes 

-24e-'^/^(£F)^f + 2a(^F)mf - log ^-^^ 



^7 



(3.4) 



where we denoted b = aj'^ to facilitate comparison with the "standard" PNJL model 
||lO| , pO|] . The weights of the adjoint representation are —2,0,2 and the group average in 
the quark sector is also easily evaluated. The result is most conveniently written in terms 
of the expectation value of the adjoint Polyakov loop. 



n„ 



a' + A' 



E 



d^fc 



3ii;^ + 2Tlog(l + e-^^/^) + 



VNf AG ^ J (27r)3 

+ 2riog(l - e-^y^ + e-2^^/^ + 3e-^^/^(£A>mf) 

This is the formula that we use for the analysis of the phase diagram. 



(3.5) 



3.1 Phase diagram 



X 

o 




The first thing that needs to be done is 
to fix the parameters of our model. There 
are altogether five of them: the coupling, 
current quark mass, and cutoff in the quark 
sector, and a, h in the gauge sector. Our 
method to estimate the NJL input param- 
eters was explained in Sec. 2A, so we sim- 



ply use the parameter set for QC2D estab- 



lished in Ref. |1C] and rescale the coupling 
according to Eq. ( 2.2[1| ). Also, we introduce 
an additional factor of two to account for the 
fact that we deal with one quark flavor only 
here. As to the gauge sector, we use the 
same physical input as in Ref. [0], that is, 
critical temperature in the pure gauge the- 
ory T^ = 270 MeV and the string tension 
(Ts = (425 MeV)^. These values were ob- 
tained from the three-color pure gauge theory using their scaling properties in the limit of a 
large number of colors, so they do not quantitatively precisely agree with those one would 
obtain directly from the two-color lattice gauge theory. However, this does not matter 
since we do not fit the parameters in the quark sector to lattice data. Within this paper, 
we merely wish to demonstrate the general trends as the number of colors or the quark 
representation are varied. 



Figure 1: Comparison of the expectation val- 
ues of the mean field a (dashed) and the funda- 
mental Polyakov loop (solid) in the naive (thin 
black lines) and Weiss (thick red lines) mean- 
field approximations to the pure gauge theory. 



Since we use a different potential for the Polyakov loop than in Ref. |10|, the param- 
eters a, 6 will actually take different values despite the same input for T^ and as- The 
deconfinement transition in the pure gauge theory is of second order with two colors, hence 
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Figure 2: Phase diagram of two-color QCD witli one flavor of adjoint quarks. Black solid line: 
deconfinement transition. Red dashed line: EEC transition. Blue dotted line: chiral crossover. The 
right panel zooms in the temperature scale so that the cusp in the deconfinement critical line is 
visible. 



we can expand the thermodynamic potential (^^) to second order in a, 



V 



(3.6) 



From here one concludes that a = Tj'logl2. The lattice spacing as, hence the parameter 
b, is then determined from the strong-coupling relation a = agCLs- The numerical values of 
all parameters are summarized in Tab. |^. 

The Weiss mean-field approximation employed here differs from the mean-field approx- 
imation used in Ref . [O , which we will henceforth refer to as "naive" for reasons explained 
in Appendix y. In the latter, the gauge sector potential can be expressed solely in terms 
of the fundamental Polyakov loop and it reads, 



j^naive 

V 



-6r[24e-'*/^4 + log(l 



(3.7) 



cf. Eq. (|3.4| ). It is therefore mandatory to compare the results obtained with the two 
approaches. We do so within the pure gauge theory. The expectation values of the fun- 
damental Polyakov loop and the mean field a are shown in Fig. |l[^ It is obvious that the 
results for the Polyakov loop are not sensitive to the particular implementation of the gauge 
sector as long as the parameters are adjusted to reproduce the same physical observables. 
Figure || shows the phase diagram of aQC2D with one quark flavor in the plane of 
temperature and quark chemical potential. The deconfinement transition associated with 
the breaking of center Z2 is denoted by the black solid line, while the BEC transition at 
which the baryon number U(1)b is broken is indicated by the red dashed line. In addition to 



*Note that there is no a in the naive mean-field approximation. The values plotted in Fig. hi were 



obtained by inverting the relation (3.3) 
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these two sharp phase transitions, there is a smooth crossover associated with the melting 
of the chiral condensate. Its position, shown in the left panel of Fig. ^ by the blue dotted 
line, is defined here by the maximum temperature gradient of a. In the chiral limit, this 
also becomes a sharp second-order phase transition. As expected, it does appear at a 
temperature much higher than that of the deconfining transition (T^ = 270 MeV, while 
T^ = 780 MeV so that T^/Td = 2.89). However, the precise value of this temperature as 
determined by our model is strongly affected by the cutoff, as is discussed in more detail 



in Sec. 4.1 
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The temperature of the deconfining 
transition depends on the chemical poten- 
tial extremely weakly, even less than in 
QC2D ll^. The reason apparently is that 
the adjoint quarks are neutral with respect 
to the center symmetry. The behavior of 
the transition lines in the vicinity of their 
"intersection" will be analyzed in detail in 
the following subsection. Finally, the BEC 
transition at zero temperature occurs at 
// = 92 MeV, which is in a good agree- 
ment with the fact that the mass of the 
pion/diquark multiplet in the vacuum is 
m^ = 184 MeV within our parameter set. 

As an illustration of the solution of the 
gap equations, we plot in Fig. |3| the con- 
densates at ^ = 100 MeV as a function of 

temperature. One can clearly see the effect of the suppression of thermal quark fluctuations 
in the confined phase: the condensates are nearly constant for T < T^. 
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Figure 3: Condensates in aQC2D a.t fi — 

100 MeV as a function of temperature. The chiral 
condensate a (black solid line), diquark condensate 
A (red dashed line), and the fundamental Polyakov 
loop (blue dotted line) are shown. 



3.2 Tetracritical point 

At the point where the two second-order transition lines cross, the system exhibits tetracrit- 



ical behavior |14|. Here we will analyze the details of the phase transitions in the vicinity 
of the tetracritical point using the GL theory. The thermodynamic potential depends on 
three mean fields, a, a, A. Only two of them, a and A, comprise order parameters for spon- 
taneous symmetry breaking of an exact symmetry (unless we consider the chiral limit) . In 
order to construct the GL free energy, one therefore needs to eliminate a in favor of a, A 
by means of its gap equation. Around the tetracritical point, we can then perform a double 
Taylor expansion of the total thermodynamic potential, 17 
and U(1)b symmetries, it depends just on the squares of the mean fields. 



n"^ + ng. Thanks to the Z2 



^}^L^ = t^a' + 6aA2 + i [A„,(a2)2 + 2X^WA^ + Aaa(A2)2] 



(3.8) 



The effective quartic couplings are determined by the second total derivatives of the ther- 
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modynamic potential, 

1 (Pn 



Vd{a 



2\2 



A 



CfA 



1 d'^n 

Vda^dA^ ' 



Aaa 



1 d^r? 

Fd(A2)2 



(3.9) 



evaluated at a = A = 0. These total derivatives are in turn given in terms of the partial 
derivatives of the thermodynamic potential as a function of all three mean fields, 



d^f) 



d^n 



d'^n fd'^VL 



dXjdxi dxidxj dxida \ da'^ J dadxj ' 



d^n 



(3.10) 



where Xi stands for a^, A^. In order to evaluate the GL quartic couplings, we need to know 
six second partial derivatives of the thermodynamic potential, 
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In order to see how the two condensates affect each other close to the tetracritical 
point, consider the general GL functional with two order parameters i?^i^2 and assume it is 
constrained to have the form 



V 



bi^l + 620^ + -{Xn4 + 2\i2(t)l4>l + A2202) • 



(3.12) 



[In our case, all other terms are prohibited by the Z2 and U(1)b symmetries.] The phase 
diagram of such a model is depicted in Fig. ^ If only one condensate were present, the 
position of the phase transition would be determined by the point where the respective h 
coefficient changes sign. However, when both condensates are present, the transition lines 
shift. This is most easily seen from the expression for the nontrivial solution to the gap 
equations with both order parameters being nonzero. 



-A22&1 + A1262 



All A 



11/^22 
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Al2fol — A1162 
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AiiA 



11^22 



A2 



(3.13) 
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We can therefore see that the size of the region with both condensates depends on the sign 
and magnitude of the offdiagonal coupling A12. 
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Figure 4: Schematic phase diagram of the Ginzburg-Landau theory with two order parameters. 
Thick hnes denote second-order phase transitions. The labels indicate which order parameters take 
nonzero values in a given phase. 



The description of the phase transitions based on the GL theory is universal and 
model independent as long as it captures the correct degrees of freedom and symmetries. 
A nontrivial task in general is to find the mapping of the (61, 62) plane displayed in Fig. ^ to 
the physical observables, in our case the temperature and chemical potential. Even though 
this is in principle possible with our PNJL model, in the present work we performed just 
a basic compatibility check. Evaluating the GL coefficients for our parameter set using 
Eq. ( pip , one finds that A„„ « 2.3 x lO'^A^, A^a » 5.7 x IQ-'^A^, and Aaa ~ 9-7 x lO'^. 
The offdiagonal coupling is positive which means that the two condensates "repel" each 
other as in the left panel of Fig. ^. However, since the GL couplings are numerically very 
small, the angles between the critical lines hardly change at the tetracritical point. The 
slight deflection of the BEG transition line is visible in the left panel of Fig. 0. That the 
same happens to the deconfinement line is made manifest by the detail of the critical line 
shown in the right panel of Fig. y. 

3.3 Casimir scaling 



The Casimir scaling hypothesis j37, ^] states that the color-singlet potential between a 
static quark and antiquark at intermediate distance is proportional to the quadratic Gasimir 
invariant, C2(7^), of the representation TZ of the quarks. This statement is exact at two- loop 
order in perturbation theory |^^ and receives corrections only at three- loop order [^^. At 
the same time, there is compelling evidence from lattice simulations that it holds to a high 
accuracy even in the nonperturbative regime ||I^, ^, ^, 43|. It may thus provide a handle 
to understand the nonperturbative behavior of QGD-like theories, and as such should be 
a necessary ingredient in any model attempting to mimic QGD (thermo) dynamics fi^ . 

In the PNJL model, one cannot directly access the confining potential feature of QGD. 
However, the scaling of the static potential implies an analogous property of the expectation 
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Figure 5: Expectation values of the Polyakov loops in various representations as a function of the 
fundamental Polyakov loop in the case of two colors. Boldface numbers indicate the "spin" j of the 
representation. Left panel: unsealed Polyakov loops. Right panel: Casimir-scaled Polyakov loops. 
For convenience, we take the C2{F)/C2(TV) power of the expectation values of the Polyakov loops 
so that the fundamental loop is left intact. 



values of the Polyakov loops [^, 45 1: the quantity {d-Ti)^ ^^ should be independent of 
the representation TZ. Since we have the analytic formula (|3.3| ) for the expectation values 
of all Polyakov loops in two-color QCD, where one has simply C2(j) = j{j + 1), we can 
easily check to what extent Casimir scaling is satisfied by our model. 

Note that the expectation values of all Polyakov loops are expressed in terms of the 
mean field a, which can in turn be traded for the fundamental loop. In Fig. |5| we therefore 
plot the expectation values of the Polyakov loops in selected representations against that in 
the fundamental representation [^, ^. Comparing the left and right panels that display 
the unsealed and scaled Polyakov loops, we can see that the Casimir scaling is very well 
reproduced as the value of the fundamental loop approaches one, which corresponds to 
high temperatures. It becomes worse at low temperatures, where the nearest-neighbor 
interaction model ( |2.1[) is too oversimplified. Lattice data that hint at almost exact scaling 
even at low temperatures can be reproduced more satisfactorily once we add more terms 
including higher representation Polyakov loops in Eq. (p.lD [19[. 



Within our model, we can check even analytically how well Casimir scaling is satisfied 
at high temperatures, and hence, at high values of a. Carrying out the Taylor expansion 
of Eq. (|3.3|) around a = -|-oo, one finds 



p \i/ia+i) 



1 1 f +j -\ .^ / 1 



a 4a-^ 12a^ \oi 

We can see that Casimir scaling is only violated at the fourth order of the expansion. 



(3.14) 
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One important observation regarding our results in Fig. ^ is that they are based just 
on the group average ( |2.8D and do not make any reference to the quark sector of the model. 
Therefore, they apply equally well to two-color QCD with quarks in any representation 
as well as to the pure gauge theory. In particular, the same curves hold even for nonzero 
chemical potential, which provides us with a unique opportunity to study Casimir scaling 
at nonzero baryon density. The quark sector will just affect the dependence of the mean 
field a on the temperature and chemical potential, and therefore the speed at which the 
curves are traversed as T and n vary. 

4. Three colors 

For three colors, the group integration is performed with the measure 

dL = ^^^ [sin(ei - 02) - sin(20i + 62) + sin(0i + 2^2)] ' ■ (4.1) 

Three-color QCD with fundamental quarks has a charge conjugation invariance, which is 
implemented in the PNJL model by a simultaneous change 6i — )• —9i, /i — )■ — //. Therefore, 
at any fixed nonzero chemical potential this charge conjugation invariance is explicitly 
broken. As a result, the expectation values (£f) and (ip) split. At the same time, the 
mean- field /3 becomes nonzero [ pO[] . 

On the other hand, the situation in aQCD is different. Thanks to the reality of 
the gauge group representation, the nonzero weights appear in pairs with opposite sign. 
Consequently, the theory is invariant under separate charge conjugation in the quark and 
gluon sectors. The charge conjugation invariance in the gauge sector guarantees that 
the Polyakov loop in a given (e.g. fundamental) representation and its complex conjugate 
always have the same expectation value. We may therefore dispense with the mean field 
P, which greatly simplifies the group integration. In the gauge sector one can still obtain 
an analytic expression for the thermodynamic potential, albeit in the form of an infinite 
series [^]. One defines a function 

+00 
F{a)= ^ det/^+i_,(a), (4.2) 

m=— 00 

where the determinant is taken with respect to the indices i,j. One then finds the following 
expression for the thermodynamic potential. 



l^y 



3 



9_^ - _6p-"/^ 
TV 



F'ia) 



F'(a) 
+ a-^ - log F{a), (4.3) 



F(a) 
and the expectation value of the fundamental Polyakov loop. 

The derivation of this formula is deferred to Appendix |C| where it will be generalized and 
used to write analytic expressions for the expectation values of all Polyakov loops. 
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a [MeV] 


61/3 [MeV] 


A [MeV] 


G [GeV-2] 


mo [MeV] 


542.1 


333.2 


651 


8.51 


5.5 



Table 2: Model parameters for three-color QCD with adjoint quarks. 

The eigenvalues of the Polyakov loop in the adjoint representation are 1 [(A'c — l)-tinies 
degenerate] and e'^^'"^-'^ for all pairs i ^ j- The logarithmic term in Eq. (p. 18) becomes 



/ ^= \ 

2\ogl{l + x)^^-^W[l + 2xcos{ei-ej) + x^]) , (4.5) 

\ i^j I mf 

where we abbreviated x = e~^ki'^ . Specifically for three colors this is equal to 

2 log| (1 + xf [l + 2xwi + x^(3 + 4a;2) + 4x3 (wi + 2^3) + x^(3 + AUJ2) + 2x^wi + x^] | • (4.6) 
Group integration reduces to evaluation of three averages, 

LOi =(cOs(6'i - 62) + COs(6'2 - 6*3) + COs(6'3 - 6'i))inf , 

UJ2 =(cos(6'i - 62) cos{6^ - Oi) + cos(6'2 - 6*3) cos(6'i - 62) + cos(6'3 - 6*1) cos(6i2 - d^))^i , 

UJ3 =(005(6*1 - 62) 005(6*2 - 6*3) 005(6*3 - 6'l))mf • 

(4.7) 

These can be performed independently of the value of x, so the evaluation of the quark 
thermodynamic potential factorizes into a one-dimensional momentum integral and a two- 
dimensional group integration. The latter can be performed either numerically or even 
analytically in a fashion similar to Eq. (|4.2|), as sketched in Appendix y. 

4.1 Phase diagram 

Again, we fix the parameters for the subsequent numerical computations first. The pa- 
rameter a is determined by the deconfinement temperature T^ in the pure gauge theory. 
With the thermodynamic potential (O), this corresponds to e"'*'^^ = 0.13427. Demand- 



ing T^ = 270 MeV, this yields a = 542.1 MeV. The parameter b is in turn obtained from 
the physical string tension ag = (425 MeV)^, as in the two-color case. In the NJL sector, 
we use the parameters of the two-flavor model with fundamental quarks, A = 651 MeV, 
G = 5.04 GeV~^, niQ = 5.5 MeV, fitted to reproduce the pion mass and decay constant 
and the chiral condensate in the vacuum (see, for instance, Ref. [^3|). The coupling is 
rescaled by the factor 27/32 in accord with Eq. (2.20), and an additional factor of two to 



account for the fact that we deal just with one flavor here. The values of all parameters 
used in our calculations are summarized in Tab. 0. 

As a basic cross-check we again evaluated first the deconfinement and chiral restoration 
temperatures (in the chiral limit) at zero chemical potential. The values T^ = 270 MeV 
and T-^ = 663 MeV yield the ratio Ty^/T^i = 2.46. This is quite far from the value ~ 8 
measured on the lattice [^, |2^. (Note that in Ref. [12| the lattice value of this ratio was 
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Figure 6: Phase diagram of three-color QCD with one flavor of adjoint quarks. Black solid line: 
deconfinement transition. Red dashed line: BEC transition. Blue dotted line: chiral crossover. The 
right panel zooms in the chemical potential and temperature scales so that the two tricritical points 
are discernible. 

achieved by tuning the parameters of the model.) However, one should keep in mind that 
we made just a rough estimate of the NJL coupling G and cutoff A, on which the chiral 
restoration temperature depends very sensitively. In principle, one could use the lattice 
value for the ratio T^/T^ as an input in the model. Nevertheless, one cannot really hope 
to describe the chiral restoration in a quantitatively satisfactory manner within our model. 
The first reason is that at such high temperatures, the calculation of the thermodynamic 
potential is plagued by cutoff artifacts. (We regulate the whole quark contribution to the 
thermodynamic potential, including its finite thermal part.) The second reason is that 
the PNJL model ceases to be physically appropriate at temperatures about two to three 



times Td p5[, since it does not capture the correct gauge degrees of freedom, that is, the 
deconfined transversely polarized gluons. We are therefore just content with demonstrating 
that QCD with adjoint quarks indeed features a large splitting of the deconfinement and 
chiral restoration temperatures. 

The phase diagram of aQCD determined within our PNJL model is shown in Fig. |6|. 
While on the large scale it looks the same as the phase diagram of aQC2D in Fig. ^, there is 
a marked difference in the topology as one zooms in the neighborhood of the "intersection" 
of the deconfinement and BEC transition lines. Since the deconfinement transition is now 
first order, the BEC critical line is broken, meeting the deconfinement line at two tricritical 
points. Thus, there is a narrow range of chemical potentials in which, as the temperature 
is increased, the diquark condensate rather unusually disappears in a first-order phase 
transition. 

4.2 Casimir scaling 

Any irreducible representation of SU(3) can be uniquely characterized by a pair of positive 
integers {p, q) that determine the highest weight of the representation in the basis of the 
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Figure 7: Expectation values of the Polyakov loops in various representations as a function of the 
fundamental Polyakov loop in the case of three colors. Boldface numbers indicate the dimension 
(and possibly the symmetry) of the representation. Left panel: unsealed Polyakov loops. Right 
panel: Casimir-scaled Polyakov loops. For convenience, we take the C2(F)/C2(7^) = l/dvj power 
of the expectation values of the Polyakov loops so that the fundamental loop is left intact. For the 
sake of clarity, the labels are not shown in the right panel. The color assignment of the lines is the 
same as in the left panel. 



fundamental weights. The triplet representation thus corresponds to (1,0) and its complex 
conjugate to (0,1). The dimension of a general irreducible representation is dim(p, g) = 
^{p+l){q+l){p + q + 2) and the value of the quadratic Casimir invariant (up to a common 
prefactor) is C2{p,q) = |(p^ + pq + q^) -\- p + q [47|. Following Refs. |2^, Q, we have 
calculated the expectation values of the Polyakov loops in the lowest few representations, 
satisfying p + g < 4. The results are shown in Fig. 0. 

As before, these results are largely independent of the quark content of the theory. 
The only assumption made is that the mean field /3 is zero so that there is a one-to-one 
correspondence between the mean field a and the expectation value of the fundamental 
Polyakov loop. Thus, the plots in Fig. apply to three-color QCD modeled by the action 
( |2.1| ) with quarks in any representation at zero chemical potential. Once the quark rep- 
resentation is (pseudo)real, the same results are valid even at nonzero chemical potential. 
As compared to the two-color case shown in Fig. |5|, the scaling violation seems to be sig- 
nificantly smaller for three colors. However, this observation is somewhat misleading since 
even the unsealed Polyakov loops show smaller depletion compared to the fundamental 
loop in the three-color case. 



5. Conclusions 

In the present paper, we worked out a description of the thermodynamics of QCD-like 
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theories at nonzero temperature and baryon chemical potential based on the PNJL model. 
To mimic the gauge sector we used a lattice spin model with nearest-neighbor interactions 
whose parameters are fixed with the help of the strong-coupling expansion of the full lattice 
gauge theory. The quark sector was modeled using the standard NJL model. 

We derived a simple mean-field expression for the thermodynamic potential, which 
is applicable to QCD-like theories with any number of colors and with quarks in any 
representation, as long as this representation is (pseudo)real. In a sequel to Ref. |11], we 



constructed the NJL Lagrangians for the two classes of QCD-like theories, denoted as type I 
and type II. 

We showed at hand of the example of QCD with adjoint quarks that the Weiss mean- 
field approximation to the lattice spin model used here is superior to the naive mean- 
field approximation, commonly employed in literature, which leads to a thermodynamic 
instability. The Weiss mean-field approximation also allowed us to derive the expectation 
value of the Polyakov loop in an arbitrary representation. The results are given in an 
implicit form applicable at all temperatures and chemical potentials, which enables us to 
study Casimir scaling in hot and/or dense matter. 

As a concrete example, we studied the phase diagram of QCD with adjoint quarks 
of two and three colors. We confirmed that in adjoint QCD the critical temperature for 
chiral restoration is much higher than that of deconfinement, both being well-defined phase 
transitions associated with spontaneous breaking/restoration of an exact symmetry (the 
former in the chiral limit). We checked the model-independent prediction that the phase 
diagram of aQC2D features a tetracritical point. On the contrary, in the phase diagram 
of aQCD the second-order BEC transition line is interrupted and meets the first-order 
deconfinement line at two tricritical points. 

It is worth emphasizing that while fine numerical details of our phase diagrams de- 
pend on our guess for the model parameters as well as on the particular way quarks are 
implemented, their qualitative features are largely based on symmetry and thus model- 
independent. Moreover, our results for Casimir scaling do not depend on the quark sector, 
in particular on the choice of the NJL parameters. They can therefore be understood as a 
direct test of the lattice spin model with nearest-neighbor interactions. Once a model for 
the quark sector is introduced, they give a prediction for Casimir scaling of Polyakov-loop 
expectation values in the whole phase diagram. 

In view of the recent lattice data E^, 49 1, the study of QCD-like theories offers a unique 



opportunity to gain more insight in the nature of strongly interacting dense matter. Even 
though all physical predictions eventually have to be made within the full gauge theory, 
we hope to have demonstrated that the PNJL model provides a versatile tool suitable for 
quick calculations and qualitative checks of the robust properties of the theory. 
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A. Fierz transformation of the current— current interaction 

Consider a fermionic field ip transforming in a representation TZ of the symmetry group. In 
NJL-like models, one deals with contact four-fermion interactions of the type X]a(V^-'^a V')^! 
where F^ is a set of matrices that project out a particular irreducible component A of the 
product representation 1Z®TZ. The Fierz rearrangement of the four-fermion interaction is 
equivalent to the group-theoretical identity 

a B b 

where the coefficients Cj,b depend only on the representations A,B. In order to fix the 
effective coupling in the meson channel, we need not evaluate the Fierz coefficients for all 

B. All we need to know is the coefficient for the one-dimensional representation B = I, 
which is always contained in the product TZ®TZ. 

Setting T^ = 1, the coefficient C_ai is projected out by multiplying Eq. ( |A.1| ) by SuSjk, 
which yields 

ETr(r^r^) 

Cja = ^- — -— ^ . (A.2) 

■^ (dim7^)2 ^ ^ 

In particular for A = X this leads to Cxi = l/dimT^. This explains the ^/Nf factor in 
the effective NJL couplings derived from the current-current interaction ( p.l9|) : both the 
original interaction as well as the term {'tptp)'^ whose coefficient we calculate are in the 
flavor-singlet channel. Likewise, the Fierz transformation from the Lorentz-vector channel 
to the Lorentz-scalar channel has the Fierz coefficient one. 

The color structure of the current-current interaction (2.19) is such that A corresponds 



to the adjoint representation, that is, F^ = Tan are the generators of the color group in 
the representation TZ of the quark fields. The Fierz coefficient ( |A.2| ) then reduces to C_4i = 
6*2(7^)/ dim 7^. Specifically for the SU(A^) group, once the generators in the fundamental 
representation are normalized as TriTaFThp) = ^6abi one finds C2(F) = (A^^ — 1)/(2A^) for 
the fundamental and C2(A) = A^ for the adjoint representation |^]. This concludes the 



derivation of the effective NJL couplings Gp and Ga given below Eq. ( 2.1 



B. Gauge group averaging with continuum quarks 

In this appendix we justify our prescription ( |2.18| ) for adding quarks to the lattice model 
of the gauge sector. In contrast to Eq. ( 2.1^ ), the authors of Ref. [^] calculated the quark 
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thermodynamic potential O^ in the mean-field NJL model with a constant background 
gauge field and set {^q)mi as the quark contribution to the thermodynamic potential. 

To start, let us emphasize that any attempt at adding continuum quarks to a lattice 
gauge model is at best heuristic. For a proper treatment one would need to discretize the 
quark action as well, thereby losing the computational simplicity of the mean-field NJL 
model. With this in mind, below we provide a qualitative argument why Eq. ( |2.18| ) is a 
reasonable approximation. 



Imagine adding quarks to the lattice model (2J); the full action then formally reads 
S^ = ^g + '4)'3>il)^ where Q} is the Dirac operator including the background gauge field the 
quarks interact with. The full partition function of the system is obtained as 

J^= f dLdi^di^e--^ = f dLe-^^det^. (B.l) 

Using the same trick of introducing the Weiss mean-field action as in Sec. |2|, this leads to 

2^ = U-^'^'9-'S^'--^i) det q\ f dL e--^-f . (B.2) 

This expression is still exact and includes all correlations between the gauge and the quark 
sectors. However, to evaluate it numerically would be very demanding. We therefore 
perform a mean-field approximation by setting 

L-i-^a-^mi) det S)\ ^ e-^-^s-'-^-f^-f (det &)^{ . (B.3) 

This is equivalent to the Weiss mean-field approximation introduced in Sec. |2| plus 
neglecting the correlations between the gauge and quark sectors.^ The full thermodynamic 
potential is then given by the gauge part (^) augmented with — Tlog(det ^)mf- One can 
therefore see that averaging the determinant of the Dirac operator is more natural than 



averaging its logarithm. However, Eq. ( 2.18 ) commits one more approximation: it neglects 
correlations between modes of different momentum and spin. While the former is naturally 
incorporated in Eq. ( p. 18 ) by the momentum integral, the latter has to be imposed by 



hand (by adding the power 1/2 to the argument of the logarithm) in presence of a diquark 
condensate, since this ties together quarks of opposite spin. Somewhat ambiguous as this 



procedure is, it does reproduce the prescription of Abuki and Fukushima [2C] when A = 0, 
and, unlike other prescriptions, it leads to a thermodynamically consistent potential il.q as 
will now be discussed. 

Let us start rather generally by addressing the following question: why have we used the 
complicated-looking Weiss mean-field approximation instead of the simple "naive" one?^ 
To find the answer it is useful to understand the relation between the two approximations. 



Let us write the Haar measure (2^) as 

Nc-l 



dL = H{e) W dOi . (B.4) 



=1 



^The lack of correlations, in particular the feedback from the dense quark matter into the gauge sector, 
makes the usual PNJL model rather trivial in the region of cold dense matter. It would be interesting to 
see to what extent these correlations can be taken into account within the present model. 

®We are indebted to Kenji Fukushima for clarifying this point at the initial stage of the project. 
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The group integral of a given function f{0), weighted by the mean- field action, can then 
be expressed as 

/ dL/(0)e-^'"f = / JJ d0i/(^)e--^-f+'°e^(^). (B.5) 

1=1 

While in the Weiss mean-field approximation this group integral is evaluated exactly, the 
naive mean-field approximation can be obtained by picking the contribution of the saddle 
point of the "action" y^i — log H{9). Indeed, let the saddle point, depending on a,/3, be 
^mf- Then the above integral is approximated by J(0j^f)e~ ™'^ """"*" ^^^™". The average 
of any function of the Polyakov loop is thus simply 

(/(^))mf = /(M • (B.6) 



Then, in the gauge thermodynamic potential (|2.9| ), the Weiss mean fields a, /3 drop out 
and the result depends only on Omf, 

Q naive 3 

'r^y " = -2{d - l)N^e-''/^£F{e^f)iUOmi) - log//(M • (B.7) 

In some particular cases, it can even be expressed solely in terms of the traced Polyakov 
loop. 

Let us now for simplicity assume that the chemical potential is zero so that there is no 
pairing and the Polyakov loop and its complex conjugate give rise to the same expectation 
values. The quasiparticle contribution to the quark thermodynamic potential ( 2.1^ ) with 
quarks in the representation TZ of the gauge group then reads 



/(W [^fcd™^ + ^Trlog(l + L7^e-^'=/^)] . (B. 



The fundamental and adjoint Polyakov loops are related by Tr^A = jTrLpI — 1, hence 
the same relation holds for their expectation values in the naive mean-field approximation. 
This means that at low temperature when the fundamental Polyakov loop goes to zero, 
the adjoint loop should become negative. Disregarding the obvious disagreement of this 
conclusion with lattice simulations, it would moreover be a disaster for the mean-field PNJL 
model. Indeed, at low temperatures, 

Tr log(l + Lne-'"/^) ^ e''"/^ Tr Ln . (B.9) 

A negative value of the Polyakov loop would thus imply that the quasiquarks would give a 
negative contribution to the pressure, leading to a thermodynamic instability. We conclude 
that the naive mean-field approximation cannot be applied to QCD with adjoint quarks. 

We will now show that a similar, albeit milder, instability occurs when one defines 
the quark contribution to the thermodynamic potential by taking (ilg)mf- For the sake of 
simplicity we focus on aQC2D at low temperature. The mean field a is then strictly zero 
(deconfinement is a sharp phase transition for adjoint quarks) and the average of the quark 
thermodynamic potential is easily evaluated using the integrals (14) of Ref. |18|. In accord 
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with the general expression ( |2.18| ) (with swapped logarithm and averaging operations), one 
finds 

2 (log[(l + x)(l + 2xcos 20 + x^)] )^^ = 2[log(l + x) - x] , (B.IO) 

where x = e~ fc' . Even though the leading term, linear in x and proportional to (Tr LA)mf) 
now vanishes, the total quasiquark pressure is still negative. This negative contribution is 
numerically small, yet it makes the thermodynamics in principle ill-defined. 

It is easy to see that this problem does not arise when the group average is taken inside 
the logarithm as in Eq. ( |2.18| ). Then at low temperature when a = 0, one gets instead of 
Eq. (CT) 

2 log((l + x)(l + 2x cos 20 + x^))^f = 2 log [(1 + x){l - x + x^)] = 2 log(l + x^) . (B.ll) 

The pressure is now strictly positive and even looks like a pressure of noninteracting 
fermionic quasiparticles with energy 'iE^.. 

One comment is appropriate regarding the last claim. In the PNJL model for physical, 
three-color QCD with fundamental quarks, one observes the same behavior at low temper- 
ature. More precisely, the mean field a is never strictly zero at any nonzero temperature, 
so the quark contribution to the pressure is proportional to log(l + 3x£f + 3x^£p + x^). At 
low temperature when the Polyakov loop is suppressed this reduces to log(l + x^), which 
is usually interpreted as a manifestation of the fact that one needs three quarks to create a 
color-singlet state. This observation suggests that the PNJL model is a natural framework 
for a description of the quarkyonic phase in cold dense quark matter [^, 52, ^|. However, 



as Eq. ( |B.llD clearly shows, this is somewhat misleading: the same low-temperature be- 
havior of the pressure arises in two-color QCD with adjoint quarks, so it does not directly 
reflect the number of quarks needed to construct a color singlet. 

A second attempt at interpreting log(l + x^) might be that both examples of three- 
color fundamental and two-color adjoint quarks are governed by the dimension of the 
representation. However, in two-color QCD it is easy to calculate the same quantity with 
quarks in higher representations, showing that there is no simple general relation between 
the representation and the form of the low-temperature pressure. For instance, in aQCD 
below the deconfinement temperature, the coefficients ti;i^2,3 take on the values oji = —1, 
u}2 = 0, wa = 1/8. Consequently, the quark pressure is proportional to 2 log(l+x^+x^+x®). 

C. Group integration for SU(A^) 

In this appendix we show that some of the group integrals can be performed for arbitrary 
A^ [54, l5l]. (For the sake of legibility, we abbreviate Nc as N.) Let us define the generating 



function 

/ N _ > 

^(z,z) = me^'="''e^""'*M . (C.l) 



In order to calculate it, we write the mean-field action ( p^q ) for one lattice site as 

N 

y^i=-J2{acosei + il3smei). (C.2) 



i=l 
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Furthermore, we use the fact that the Haar measure ( |2.6D may be written as a square of a 
Vandermonde determinant, 



N 



dL = l[ dOi 5{ei + ■■■ + dN)e,,...,^e,,...jJ'^^'^~^''^ ■ ■ ■ e^'^^^^-i^) . (C.3) 

The last trick is to express the (periodic) 5-function in terms of its Fourier series, 

-, +00 
5{6i + --- + 6n) = — Y e^™(^i+-+^^) . (C.4) 

ATT f ^ 
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m=— 00 



The integration over the angles Oi now completely factorizes in terms of a single master 
integral, 

2vr Jo 

For real u and pure imaginary v, v = iw, which is the case if /3 = 0, the master integral 
can again be expressed with the help of the modified Bessel function, 

Xiu,iw) = ( ;!Z^ y/„(x/n^ + i.2) . (C.6) 

The final formula for the generating function ( p.lD reads 



y^ det 3^m+i-j {a + z + z,(3 + z - z) 

^(^,^-) = — ^^ . (C.7) 

y^ det. 9^+1- j (a, P) 



m=—oo 



Looking back at Eq. ( C.l ) one sees that expanding the exponentials, the Taylor co- 
efficient of the z^z" term resums all eigenvalues of the Polyakov loop in the F*" (8) F" 
representation, F being the fundamental one. That is, one has 






(C.8) 

z=0 

2=0 



The expectation values of Polyakov loops in all irreducible representations can be obtained 
from this formula by simply observing that the (traced) Polyakov loop in a direct sum of 
two representations is equal to the sum of the loops in these representations. 

Let us remark here that the thermodynamic potential of the three-color pure gauge 
theory ([4. 3D can be derived using the same argument, and the group integrals involved 



are special cases of those considered above. Indeed, the function F{a) ( [4. 2]) equals the 



denominator in Eq. ( C.7 ) at /3 = up to a trivial numerical prefactor. Changing this 
prefactor just shifts the thermodynamic potential by a constant, and noting that F{0) = 1, 
it can be fixed by demanding that fi^ = for q = 0. 
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A more compact formula can again be obtained for the special case of two colors. 

The one-dimensional group 



ze'" 



Then, we can set /3 = and z = 0. Also, ni=i ^ 
integration can be performed directly and one finds 



„2z cos t 



^{z) 



Ioi2a + 2z) - I2i2a + 2z) a Ji (2a + 2z) 



Jo (2a) -l2 (2a) 



a + z /i(2a) 



(C.9) 



While the latter expression is more compact, the former is more convenient for taking the 
derivatives in order to extract the expectation values of the Polyakov loops. 

Finally, let us show that even the averages ( [4.7D can be expressed analytically in terms 



of a series of modified Bessel functions [56], and thus speed up the numerical evaluation of 
the thermodynamic potential. Using trigonometric identities, these averages can be written 
as a linear combination of terms of the type 



Kabcia) = (e'( 



i+b92+c93)'. 



mf ) 



(CIO) 



where a, b, c are integers. Using the same trick of rewriting the Haar measure as a Van- 



dermonde determinant and introducing the periodic (5-function as in Eqs. (|C.3| ) and ( C.4 ) 
this becomes 



Kabc{a) 



1 



6F(a) 



Yl Yl ^vk 

n=—ao ij,k=l 



Im+i-l+a{a) Im+i-2+a{oi) /m+j_^3+a(a) 
Im+j-l+b{oi) Im+j-2+b{oi) Im+j-3+b{oi) 
Im+k-l+cia) Im+k-2+c{oi) /m+fc-3+c(a) 



(C.ll) 
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